STEADY WATER WAVES WITH MULTIPLE CRITICAL LAYERS: 

INTERIOR DYNAMICS 
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Abstract. We study small-amplitude steady water waves with multi- 
' ' ^ pie critical layers. Those are rotational two-dimensional gravity-waves 

I . propagating over a perfect fluid of flnite depth. It is found that arbitrar- 

' ily many critical layers with cat's-eye vortices are possible, with different 

^ structure at different levels within the fluid. The corresponding vorticity 

, ^ ■ depends linearly on the stream function. 

Oh! 

<: 

I 1. Introduction 

■ This paper is concerned with steady, periodic, two-dimensional gravity- 
II I waves of permanent shape and velocity. Famous among these are the Stokes 

1^ ' waves; symmetric waves with a surface profile which rises and falls once in 

. every minimal period [2, 30, 31, 38, 41]. Of particular interest are waves with 

^ I vorticity. Vorticity plays a crucial role for wave-current interactions and in 

the formation of wind-generated waves [17, 36, 39]. The first mathematical 
, construction of a rotational free-surface fluid flow is due to Gerstner in the 

J> ' beginning of 19th century [4, 24], but it was first with the pioneering inves- 

■ tigation [12] that a modern theory of both large and small rotational waves 
I was established. This theory was extended to deep water waves in [26], to 
' small-amplitude waves with surface tension in [44], and to large-amplitude 
. waves with surface tension and stratification in [46]. 

I Many properties inherent in irrotational periodic gravity waves, such as 

■ the symmetry of the surface profile [34, 42], the analyticity of the stream- 
. . , lines [11], and the Stokes conjecture [1, 35], carry over to rotational water 

' waves [7, 11, 8, 9, 43]. A notable exception from this rule arises when one 

^ . examines fiows with internal stagnation, i.e. points where the velocity of a 

^ I fluid particle coincides with that of the wave itself. Even when the vorticity 

■ - - ' is only constant, critical layers with cat's-eye vortices arise [20]. Those are 

horizontal layers of closed streamlines separating the fluid into two disjoint 
regions, a behaviour that is not possible for irrotational waves [5]. Recently, 
the existence of such waves with one critical layer as solutions of the full Eu- 
ler equations was established [45] (see also [15] and, for a study of stagnation 
points in rotational flows, [13]). This constituted an important connection 
between the mathematical research on exact Stokes waves and the study 
of waves with critical layers within the wider field of fluid dynamics (cf. 
[3, 29, 33, 37, 40]). 
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In the very recent investigation [18] the theory from [45] was extended 
to the case of an affine vorticity distribution, yielding the existence of ex- 
act small-amphtude gravity waves with arbitrary many critical layers. It 
is our aim to give a quahtative, and to some extent also a quantitative, 
description of those waves ([18] also contains a proof for the existence of 
exact bichromatic waves with vorticity, but those are not considered here). 
This is a natural continuation of a line of recent research on the qualita- 
tive features of various free-surface flows (see, e.g., ]21], and the references 
[6, 10, 14, 16, 22, 23, 25, 27, 28] from that survey). In particular, the results 
here obtained could be used to give a description of the particle trajectories 
within linear waves with critical layers. 

The focus of this paper is, however, slightly different. Since the construc- 
tion in [18] is by bifurcation from laminar flows, for small waves it is possible 
to investigate the flow and give exact estimates on the error (see Proposi- 
tion 3.2). It turns out that waves with an affine vorticity distribution can 
be naturally divided into four wave classes, and we give the velocity fields 
and the bifurcation relations in the different cases. Theorem 5.2 provides 
the qualitative description of the wave class with multiple critical layers. 
Although there is a rich variety of flow configurations we discern two main 
scenarios (which, essentially, can be seen in Figure 2). An interesting feature 
is that, apart form the region closest to the free surface, the fluid motion 
takes place in vertically disconnected, completely flat regions (which is not 
the case in waves with a single critical layer). In Theorem 4.1 we are also able 
to give a quantitative description of the levels at which stagnation points, and 
therefore critical layers, can arise. Those results are graphically captured in 
Figure 1. 

The disposition is as follows. Section 2 describes the governing equations, 
with Section 3 narrowing in on laminar flows and the first-order perturbations 
thereof. In Section 4 we describe the four wave classes, and detail at which 
levels stagnation can occur in the different cases. Finally, Section 5 presents 
the main structure of the interesting wave class with multiple critical layers, 
and some numerical examples are given. For a quick glance at the waves, 
see the last section. 



Let (x, y) be Cartesian position coordinates, and (u, v) = {x, y) the corre- 
sponding velocity field. Here 



are 27r-periodic in the x-variable and the vertical coordinate y ranges from 
the flat bed at y = to the (normalized) free water surface at y = l + r/(t, x). 
Let p := p{t, x, y) denote the pressure, and g the gravitational constant of 
acceleration. In the mathematical theory of steady waves it is common and 
physically realistic to consider water as inviscid and of constant density [30, 
31]. The Euler equations 



2. Preliminaries 



u := u{t,x,y) 



V := v{t,x,y) 



Ut + UUx + VUy = -px, 

vt + uvx + vvy = -py - g, 
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then model the motion within the fluid. The equations 

Ux + Vy = and Vx — Uy = uj (2-lb) 

additionally describes incompressibility and the vorticity u, respectively.^ 
At the surface the conditions 

p = Po and V = rjt + ur]x (2.1c) 

separate the air from the water, po being the atmospheric pressure. Note 
that the second condition in (2.1c) states that y{t) — r]{t,x{t)) is constant 
over time, so that the same particles constitute the interface at all times. 
Similarly, no water penetrates the flat bed, whence we have 

v = at y = 0. (2. Id) 

The equations (2.1) govern the motion of two-dimensional gravitational wa- 
ter waves on finite depth. 

An important class of waves are travelling waves, propagating with con- 
stant shape and speed. Mathematically, such waves are solutions of (2.1) 
with an (x — ct)-dependence, where c > is the constant wavespeed, and we 
have restricted attention to waves travelling rightward with respect to the 
fixed Cartesian frame. Since Dt{x — ct) = u — c, it is natural to introduce 
steady variables, 

X := X — ct, U := u — c. 

We shall also write Y for y and V for v to indicate when we are in the 
travelling frame. In the steady variables the fluid occupies 

n^:= {{X,Y) eM.'^-.O <Y <l+r]{X)} . (2.2) 

Define the relative pressure P through 

p=:po + gil + P-Y). 

Since the term —gY measures the hydrostatic pressure distribution, the rel- 
ative pressure is a measure of the pressure perturbation induced by a passing 
wave. Altogether we obtain the governing equations 

UUx + VUy = -gPx, 

UVx + VVY = -gPY, , , 

Vx-Uy = W, 
with boundary conditions 

P = r?, 

y = Uvx, 
and 

V = 0, on y = 0. (2.3c) 

The problem of finding (U, V, P, rf) such that (2.3) is satisfied is known as 
the steady water-wave problem. Since is an a priori unknown, (2.3) is a 
free-boundary problem. 



on Y = l + r]{X), (2.3b) 



-'^This sign convention for the vorticity is consistent with [18]. The reader should be 
advised that the vorticity may also appear with the opposite sign in the literature. 
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The ct-problem. When r] € C^(]R), and u,v £ C'^{Qri), one can use the fact 
that the velocity field is divergence-free (cf. (2.1b)) to introduce a stream 
function tp G C^{Q) with 

■= -V and il^Y := U. (2.4) 
Define the Poisson bracket {/, 5} := fx BY — fvgx- 

Proposition 2.1 (Stream-function formulation). The water-wave problem (2.3) 
is equivalent to that 



(2.5) 



for some constants rriQ, mi, and C . 

Proof. Identify if) with U and V through (2.4). Given the regularity assump- 
tions and that f]^ is simply connected, we see that Ux + Vy = is equivalent 
to the existence of ip. The relations P = 77 in (2.3b) and y = in (2.3c) mean 
that "0 is constant on the surface and on the flat bed, just as Vx — Uy = —oo 
means that l!^tj) = —lo. It remains to show how the equations of motion 
relate to the Bernoulli surface condition and the bracket condition. 

Given (2.3) one can eliminate the relative pressure by taking the curl of 
the Euler equations. That yields 

UAV - VAU = 0. (2.6a) 

Moreover, by differentiating the relation P = r] along the surface, and using 
(2.3a), we find that 

U'^ + V'^ + 2gY = C, Y = l + r]{X). (2.6b) 

Hence (2.5) holds. Contrariwise, if {U,V) fulfil (2.6a) and (2.6b), one can 
define P up to a constant through (2.3a), and, with the right choice of 
constant, P satisfies (2.3b). □ 

Consider now the case when ipy may vanish, but Ail^y/i^Y can be extended 
to a continuous function on O^, i.e. 

a := ^ € C\n,) (2.7) 

One can then exchange the bracket condition {ip, Atp} = in (2.5) for 

(A-a)VV' = 0. (2.8) 

When a is a constant there exists an affine vorticity function 7 with 7' = —a, 
meaning that 

A^p = -'y{^p) = aip + (3, /3 G M. 
Observe that this does not rule out the existence of stagnation points Vip = 0. 
Without loss of generality we may take /? to be zero; changing it corresponds 
to changing mo and mi. The choice a = models constant vorticity and 
was investigated in [20, 45]. The next natural step is a constant but non- 
vanishing a. That is the setting of this investigation. 
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3. Laminar flows and their first-order perturbations 

Laminar flows are solutions of the steady water-wave problem (2.3) with 
r]{X) = 0. Those are the running streams for which 

U{X, Y) = Uo{Y) and V = P = rj = 0. 

We shall require that Uq € C^([0, The function Uq is the (rotational) 
background current, upon which we will impose a small disturbance: the 
system (2.3) will be linearized at a point {U,V, P,r]) = (Uq, 0,0,0), and the 
solutions of the constructed linear problem analyzed. We thus assume that 
U, V, P and i] allow for expansions of the form 

f = fo + efi + 0{e^) , as e ^ 0. (3.1) 

Here Uq is a background current as described above, and 

Vo = Po = m = 0. 

By inserting these expansions into (2.3), and retaining only first-order terms 
in e, we obtain the linearized system 

dxUi + dyVi = 0, 

Uo dxUi + Vi dyUo = -dxPi, in R x (0, 1) (3.2a) 

UodxVi = -dYPi, 

with boundary conditions 

Vi = Uodxm, 

on y = 1, 3.2b 

as well as 

Vi = on Y = 0. (3.2c) 
The following result allows us to eliminate the relative pressure from (3.2). 

Proposition 3.1. Let the background current Uq be given. Under the condition 
that 



/7r PTT 
7]i{X)dX = and / Ui{X,Y)dX = 0, ^€[0,1], 
-TT J —TV 



(3.3) 



the solutions {Ui,Vi,Pi,r]i) of (3.2) are in one-to-one- correspondence with 
the solutions Vi of 

UoAVi = U[;Vu 0<Y<1, 
(1 + UoU[,)Vi = dyVi, Y = l, (3.4) 

Vi=0, Y = 0. 

Proof. Taking the curl of the linearized Euler equations, and differentiating 
p = 1] along the linearized surface Y = 1 yields (3.4). If [Ui,Vi,Pi,rji) is a 
solution of (3.2) then Vi fulfills (3.4), and if Vi is a solution of (3.4), then one 
can find {Ui, Pi,r]i) such that (3.2) holds. One defines Ui through the first 
equation in (3.2a), and then Pi through the two last equations in (3.2a). The 
linear surface 7]i can be determined by (3.2b), and the boundary condition at 
y = 1 in (3.4) guarantees that (3.2b) is consistent with (3.2a). Notice, how- 
ever, that for a given Vi, a solution C/i is only determined modulo functions 
f{y), and rji up to a constant. We shall therefore require that the periodic 
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mean of the first-order solution equals that of the running stream, meaning 
that (3.3) holds. In particular, this implies that the solution {Ui,Vi, Pi,r]i) 
of (3.2) is unique with respect to the solution Vi of (3.4). □ 

Laminar vorticity. Now, suppose that Uq /Uq can be extended to a continu- 
ous function on [0, 1], and introduce the laminar vorticity 

ao:=^€C([0,l],M). (3.5) 

Let (cf. (3.4)) 

/ii := 1 + f/o(l)C/o(l) and ^2 := U^il). 
According to Proposition 3.1 we may then consider the system 
AVi = aoVi, < y < 1, 

fiiVi= fi2dYVi, Y = l, (3.6) 

Vi=0, Y = 0. 

In our case q = qq S M, and constant vorticity is captured by qq = 0. 

Relation to exact nonlinear solutions. In what comes we will find and inves- 
tigate four solution classes of (3.6) and thus of (3.2). Any exact solution 
of the steady water-wave problem (2.5) with Aip = aip that allows for an 
expansion as in (3.1) and adheres to the normalization (3.3) will satisfy the 
velocity fields here investigated up to an error of order in the appropriate 
space. A particularly case (Wave class 1 on page 7) corresponds to a class 
of solutions found in [18] by linearizing around a running stream with back- 
ground current Uo{Y) = asin{9Q{Y — 1) + A). Here 6q = -^/jooj and A G R 
is a parameter. Those solutions do not necessarily satisfy the normalization 
(3.3); while J^^ rji dX = the strength of the first-order background current 
may change with e. This is the reason why a depends on e in the following 
proposition, which is a consequence of the results from [18]. In accordance 
with (2.2) we let, for any small and positive constant 6, the set Q-s denote 
the part of the fluid domain where < Y < 1 — 6. 

Proposition 3.2. Let e ^ (^,f?) G (7^(0,^) x C^(M) be a solution curve found 
in [18] by bifurcation from a one- dimensional kernel of minimal period 27r. 
Pick < (5 << 1. For any e small enough, there exists a = a(e) such that 
the velocity field {U^V) = {ipY , —^^x) coincides with that of wave class 1 with 
Uq{Y) = a{e) sin(6'o(y - 1) + A) up to addition of terms 0{e^) in C'^(yL_s). 
The map e 1— t- a(e) is smooth and a(0) fulfils the bifurcation condition (4.3). 



4. Wave classes 

Even when uq is a constant, the linear system (3.6) contains a rich variety 
of solutions, including asymmetric ones (cf. [19]). We shall see that restrict- 
ing attention to the first Fourier mode of Vi still produces a wide range of 
linear waves. We thus search for a solution of the form 



= sin(X)/(y), /€C2([0,1],R). 



(4.1) 



CRITICAL-LAYER WATER WAVES 



7 



The ansatz (4.1) reduces the system (3.6) to a (trivial) Sturm-Liouvihe 
problem: 

-/" + (ao + 1)/ = 0, 

w/(l)-^2/'(l) = 0, (4.2) 
/(O) = 0, 

with ckq € M, //f + > 0, and H2 > 0. Since the case ao = has already 
been treated in [20] we restrict our attention to ao ^ 0, corresponding to 
non-constant vorticity. Define 



^0 := vlaol and 9i := y/lao + 1|. 

Using the Sturm-Liouville problem (4.2) to determine V, and then U via 
Proposition 3.1, we find that the solutions belong to one of the following 
four classes: 

Wave class 1 (Laminar vorticity oq < —1)- The solutions of (4.2) are gener- 
ated by f{Y) = sin(6liy) with 

Uo{Y) = asm{0o{Y -l) + X), 

= sin^ (A) {Oi cot(0i) - cot(A)) , (4.3) 

AG(arccot (^^) , ^) • 

Up to the first order in e, 

U{X, Y) = Uo{Y) + eOi cos X cos{eiY), 

V{X, Y)=e sin X sm{9iY) . ^^'^'^ 

Wave class 2 (Laminar vorticity ao = — !)• The solutions of (4.2) are gener- 
ated by f{Y) = Y with 

?7o(y) = asin(y-l + A), 

= sin2(A) (1 - cot(A)) , AG(f,7r). (4.5) 

Up to the first order in e, 

U{X,Y) = Uo{Y) + ecos{X), 

V{X,Y) = £Ysm{X). ^ 

Wave class 3 (Laminar vorticity —1 < ao < 0). The solutions of (4.2) are 
generated by f{Y) = sinh(0iy) with 

Uo{Y) = asm{eo{Y-l) + X), 

a"2 = sin^ (A) (6*1 coth(6'i) - ^o cot(A)) , (4.7) 

AE(arccot(^l^),.). 

Up to the first order in e, 

U{X,Y) = Uo{Y) + £6i cosXcosh(eiy), 

V{X,Y) = esinXsinh(6liy). ^^'^^ 



8 EHRNSTROM, ESCHER, AND VILLARI 

Wave class 4 (Laminar vorticity oq > 0). The solutions of (4.2) are generated 
by f{Y) = sinh(6liy) with 

Uo{Y) = asmh{eo{Y - 1)) + Xcosh{9o{Y - 1)), 

A20icoth(0i) - 1 , , , 
a= ' , (4.9) 

Up to the first order in e, 

U{X,Y) = Uo{Y) + e6i cosXcosh(0iy), 

V{X,Y) = esinXsinh((9iy). ^^'^^^ 



Stagnation. The explicit solutions allow us to determine the possible levels 
of stagnation. From the following result one obtains Figure 1. 

Theorem 4.1 (Stagnation). The following hold for the background current Uq 
for the wave classes 1-4- 
W 1. For any Yq € [0, 1) there exists ckq < —1 and A such that Uo{Yq) = 0, 
and the number of zeros of Uq in [0, 1] can be made arbitrarily large 
with the appropriate choice of a^. 

In the interval —1 — < oq < —1 the background current Uq has 

one zero at Yq = 1 — X/Oq for arccot(^?i cot(0i)/0o) < A < 

none for X G (^Oi '^)- 
W 2. C/q has one zero atY^ = 1 — X for X € (vr/4, 1] and none for X G (1, tt). 
W 3. Uq has one zero atY^ = 1 — X/Oq for arccot(0i coth{9i) / Oq) < X < Oq 

and none for X € {0q,tt). 

W 4. Uq has one zero at Yq = 1 — 9^^ arctanh (^ ;^ag^ coth(ei)-i ) f^''" "^^ — 
(^1 coth(0i) — 9q coth(0o))~^ o,'>^d none for other X ^ 0. 



Proof. The analysis is carried out separately for each wave class. 
W 1. The function 

9i cot (01 ) 



spans the real numbers (it blows up at oq = — 1 — n^vr^, n E N). We thus see 
from (4.3) that the set (arccot(0i cot{9i)/9o), vr) may be empty. But for any 
e > and n € N, there exists 6 > such that if < —1 — Ji^vr^ < oq + 6, 
then (4.3) is solvable for all A G (ejVr). The number of zeros is at least as 
large as [^o/"^] — >• co as oq — > — oo. The last proposition then follows by 
checking that Uq can have at most one zero Yq for — 1 — 7r^<a<— 1. 
W 2. Consider (4.5). The function 

TT 

A I— > sin(A) (sin(A) — cos(A)) > exactly when — < A < vr, 

and since it is bounded, the amplitude a is bounded away from 0. There 
thus exist A and Yq € [0, 1] such that Uq{Yq) = if and only iiYQ = l — X£ 

[0,1-f]. 

W 3. The right-hand side of (4.7) is positive when 
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Figure 1 . The possible zeros Yq of the background current Uq as 
a (multivalued) function of the laminar vorticity ao; those are the 
levels of the critical layers in the limit e ^> 0. For each pair (ao, Yq) 
in the shaded region, and only for those, there exists A such that 
the bifurcation condition is fulfilled and Uq has precisely one zero 
at Yq. As ao oo we have maxYo ~^ 1, but for any given ao the 
stagnation points are bounded away from the surface as e — > 0. In 
contrast, stagnation at the flat bed is possible whenever ao > — tt^. 
For ao < — 1 — tt^ more zeros appear and the situation is not as 
transparent. 



To have Yq € [0,1] with Uo{Yo) = necessarily A G [0, ^o]- The assertion 
follows from that is strictly larger than the lower bound in (4.11). 

W 4. The background current has at most one zero, and to see what zeros 
there are in [0, 1] we consider 

tanh(0o(^ - 1)) = -- = r, ^° , 

^ ^ ^' a I- A26iicoth(6ii) 

For Yq € [0, 1) the left-hand side is negative, so that we must at least have 
> (6*1 coth(6'i))~^, and a closer look yields that A^ > (0icoth(6'i) - 
00 coth(0o))~^ is required to match y > 0. The right-hand side is then 
an increasing function of A^, and 

- tanh(0o) < \^.,n^ < 



1 - A26'i coth(6'i) - 6'icoth(6ii)' 

Since also arctanh is an increasing function, the question reduces to whether 

1 / 9 

0<Yo<l-— arctanh 

9o V^icoth(6li) 

The right-hand side is positive, strictly increasing in ao, and tends to 1 as 
ao — >■ oo. □ 

5. HAMILTONIAN FORMULATION AND PHASE-PORTRAIT ANALYSIS 

In the analysis to come the region of interest is 

< y < l±0{e) cos{X), 

where the sign indicates that for each of the wave classes 1-4 one finds that 
X = may be either a crest or a trough, depending on the signs of a and A. 
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Figure 2. Two scenarios for wave class 1. Blue lines are stream- 
lines, dotted lines isoclines [oo-isocline green, 0-isocline red], and 
fat dots critical points [centers red, saddle points green]. Left: 
an uppermost critical layer with horizontal streamlines cutting 
through it as in Thm. 5.2 ii.b), and a lower critical layer as in 
Thm. 5.2 ii.a). Right: two critical layers as in Thm. 5.2 ii.a) sep- 
arated by a horizontal streamline. Note that the rotational flow 
near the bottom evolves under a "rigid lid". 



Recall that {u, v) = (x, y). In view of that X = x — ct and U = u — c, one 
similarly obtains 

{X,Y) = {U,V). 

The paths (X{t),Y{t)) describe the particle trajectories in the steady vari- 
ables, and any such solution is entirely contained in one streamline. 

Proposition 5.1 (Hamiltonian formulation). The wave classes 1-4 admit 
a Hamiltonian 

H{X,Y) := ecos{X)G{Y)+ f Uo{s)ds, (5.1) 

with 

sin(0iy), for wave class 1, 

G{Y) := I Y, for wave class 2, (5.2) 

sinh{9iY), for wave classes 3-4- 

Classes 2-4 can be dealt with as the class a = (constant vorticity) in [20, 
45], and do not yield any new qualitative results. Indeed, their appearance 
and the analysis thereof is captured within that of the interesting class 1. 

Theorem 5.2 (Wave class 1). The following hold for small- amplitude waves 
of wave class 1 (e sufficiently small). 

i. The fluid motion is divided into vertical layers, each separated from 
the others by flat sets of streamlines {(X, y*): sm{6iY^:) = 0}. 

ii. For each Y^, with C/o(y*) = there is a smooth connected part of the 
oo-isocline passing through all points (vr/2 -|- mr,Y^:), n € Z, along 
which centers (cats-eye vortices) and saddle points alternate in one 
of the following ways: 

a) when Y^, is not a common zero ofUo and sm{9i-) centers appear 
at every other X = rnr and saddle points at every other (n-|-l)7r; 

b) when Y^, is a common zero ofUo and sin(0i-) centers appear at 
X = HTT and saddle points at tt/2 -\- mr. 



Remark 5.3. Starting with the situation in ii.b) one might fix e and a, and 
then vary the zero of the background flow. The saddle point at X = 7r/2 
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Figure 3. Numerical plots of wave class 1. Colours indicate the 
strength of the velocity field. Top left \a = —20, A = 4.39, e = 
0.05]: the case of a common zero as in Thm. 5.2 ii.b). The upper 
critical layer is magnified in the plot bottom left, where the three 
centers, the two saddle points and the 0-isocline cutting through 
them are clearly visible. The bottom right plot shows the same 
flow near the bottom with only two centers as in Thm. 5.2 ii.a). 
The situation is similar to that in Figure 2, left, with the difference 
that the oo-isocline actually crosses the fiat bed in this numerical 
plot. Right [a = -20, A = 4.60, s = 0.05]: as the zero of the 
background current shifts, the horizontal 0-isocline climbs above 
the critical points in the upper critical layer and the saddle nodes 
merge with one center each, turning the situation in Thm. 5.2 ii.b) 
to the one in Thm. 5.2 ii.a). 

then continuously and monotonically approaches the center at either X = 
or X = TT, eventually merging with it and wiping it out. 

Proof. It follows from (5.1) that the fluid motion is 27r-periodic and sym- 
metric around the vertical X = axis. It therefore suffices to investigate the 
strip <X <TT. 

i) . The velocity field is 

X = Uo{Y) + eOi cos{X) cos{eiY), 
Y = esm{X)sm{eiY), 

whence the 0-isocline, defined as the set where y = 0, is given by the vertical 
axes X = mod vr and the horizontal lines where sin(^iy) = 0. 

ii) . Let y, be a zero of C/q. Then (7r/2, y*) belongs to the oo-isoc/me {(X, y) : X = 
0}. Since Uq{Y^:) = oOq ^ 0, we have that 

DyX = Uq{Y) - eel cos(X) sin(^iy) 
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= ....r r^rr.:;. (5-3) 



is nonzero at (7r/2,y*). The implicit function theorem allows us to locally 
parametrize the oo-isocline as the graph of a smooth function Yao{X) with 
slope 

eOi sin(X)cos(giy) 
U[,{Y) - eel cos(X) sin(0iy) ' 

A continuity argument yields that, for e small enough, Yqo extends to a 2ti 
periodic function on M, strictly rising and falling between the zeros of sin(X), 
and with \Yoo -Y^\ = 0{e). Hence U'q{Y^) 7^ 0. At X = nvr, n G Z, the 
graph of Yoo intersects the 0-isocline. At those critical points, the Hessian of 
the Hamiltonian is given by 



J}^H[n^,Y^\x-. 



-l)'^+iesin(0iyoo) 

{-ir+^eel smiOiY^) + U[,iY^) 



(5.4) 

X=n7z 



There are now two cases. 



a) When sin(0i-) and Uq have no common zero. From |yoo — Y*\ = C^i^) we 
find that sin(0iyoo) is non-vanishing and thus of constant sign. For e small 
enough, the Hessian (5.4) thus has one negative and one positive eigenvalue 
at every other X = rnr, and two of the same sign at every other X = {n + l)7r 
in between. The assertion ii.a) then follows from the Morse lemma [32]. 

b) When sin(0i-) and Uq have a common zero. In this case there are ad- 
ditional critical points a,t X = 7r/2 -|- nvr, n (z 1^, all similar to the one at 
X = 7t/2. There 



I)^H{tt/2,Y; 



-£6*1 cos(6lin) 

-e^i cos(0iy,) c/^(y,) 



with one strictly positive and one strictly negative eigenvalue. Hence, the 
critical point (7r/2,y*) is always a saddle point. 

We want to show that the Hessian (5.4) has two eigenvalues of the same 
sign at all critical points {mT,YoQ\x=nn)- Since the slope of yoo changes 
direction exactly at X = mr it follows that also in the case when C/o(y*) = 
sm{6iY^:) = we have sin(0iyoo|x=n7r) 7^ 0, but with 

sgnsin(6'iyoo|x=n7r) = -sgnsin(6'iyoo|x=n7r), 

all given that e is small enough. We now claim that —esm{9iYoo\x=o) a-^d 
Uq{Yoo\x=o) have the same sign (cf. (5.4)). The slope of yoo|xG(o,7r) 
determined by the sign of cos(0iyoo)/C/o(^^) same interval. We have 

cos(0i-) > when sin(^i-) is increasing, and contrariwise. The assertion now 
follows from that sin(0iy*) = . □ 
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